Analytical Solution of Cross Polarization Dynamics 
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Cross polarization (CP) dynamics, which was remained unknown for five decades, has been derived 
analytically in the zero- and double-quantum spaces. The initial polarization in the double-quantum 
space is a constant of motion under strong pulse condition QtJii + wis\ 2> while the Hamil- 

tonian in the zero-quantum space reduces to d(t)a^ under the Hartmann-Hahn match condition 
(ojii = Wis). The time dependent Hamilontian (d(t)a^) in the zero-quantum space can be ex- 
pressed by average Hamiltonians. Since[d(t')cr A , d(t")o~ A ] = 0, only zero order average Hamiltonian 
needs to be calculated, leading to an analytical solution of CP dynamics. 
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1. INTRODUCTION 

Cross polarization (CP), developed by Hartmann and 
Hahn [l[ in 1962 and later modified by Pines et al 
is one of the most important techniques in NMR. Medi- 
ated by heteronuclear dipolar interaction, polarization is 
often transferred from abundant I spins to a rare S spin 
with a polarization enhancement up to 77/75. Because 
the overall delay time is just only subject to the spin- 
lattice relaxation time of I spins, considerable NMR time 
is saved compared with a single pulse experiment. 

CP dynamics for a single crystal was first derived by 
Miillcr et al Q . For a stationary sample, the Hamiltonian 
is time independent, it can also be solved by diagonaliza- 
tion of the Hamiltonians. Under fast magic angle spin- 
ning(MAS) (cj r /2-7r > \d(t)\), the Hartman-Hahn match 
condition shift to u>u — ujis = nw r (n = ±1,2) 0] and 
the effective heteronuclear dipolar interaction becomes 
time independent [B| . Consequently the CP dynamics un- 
der fast MAS can be derived without many problems [f| , 
in particular under Lee-Goldburg (LG) [6j condition as 
demonstrated by Ladizhansky et al [7|, [8(. However, 
under conventional MAS speed, the heteronuclear dipo- 
lar interaction becomes time dependent, as many other 
quantum systems, searching for an analytic solution is 
usually not conceivable. 

So far, CP dynamics under conventional MAS speed 
is usually calculated with a spin temperature hypoth- 
esis I 9 | . Ild| . It inevitably leads to an empirical solu- 
tion [1]| . In the zero- and double-quantum spaces [l2j ]. 
which commutes with each other, the evolution of density 
matrix can be done separately in the two spaces, sim- 
plifying the calculation considerably. By means of this 
method, a number of intriguing phenomena in CP has 
been understood thoroughly. They include CP d yna mics 
of phase-shifted CP under mismatch conditions Il2l . 1 1 3j j , 
W-MOIST Q, double-quantum matched CP @,[l3, and 
adiabatic polarization transfer [l6[ • 

Up to now, CP dynamics under conventional MAS 
speed and the Hartman-Hahn match condition (u>n = 
wis) has remained unknown. In this article we show that 



this problem can be resolved in the zero- and double- 
quantum spaces together with average Hamiltionian the- 
ory [13] • The experiment and simulated results match 
well with the theoretical predictions. 

2. THEORY 

For a heteronuclear dipolar coupled IS spin pair under 
cross polarization and magic angle spinning (CPMAS), 
the Hamiltonian can be described by 



TL = Lu u I y + LJisSy + 2d(t)I z S z 



(1) 



where oju and u>is are the strengths of the I and S spin- 
locking fields, respectively, and 



d(t) = d[V2sin(2/3)cos(uj r t + 7) 
- sin 2 (/3)cos(2ui r t + 27)] 



(2) 



is a time dependent heteronuclear dipolar coupling con- 
stant with two Euler angles in the rotor-fixed frame. In 
the zero- and double-quantum spaces, the above Hamil- 
tonian can be expressed as 0, [l2Hl4j 

% A = d(t)a A + (uiu — ui\s)o- A {zero — quantum) (3) 

and 

"H s = d(t)af + (u)u + 0Jis)cr y (double — quantum), (4) 

where a A and of are the Pauli matrices. The initial 
density matrix can be expressed in terms of zero-quantum 
and double-quantum density matrices 



P (0) = Iy 



Sy 



(5) 



Because [<rf ,of] = 0(i,j = x,y,z) and [H A ,H^} = 0, 
the evolution of density matrix can be calculated in the 
zero- and double-quantum spaces separately, 



p(t) 



a A (t) + cr^(t) 
= e -iJo' HA ( t ') dt 'a A e^ « A (*')<it' 



(6) 
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Under strong pulse condition \lju + Wis] 3> \d(t)\, the 
density matrix in the double-quantum space is nearly a 
constant of motion, i.e. 



(7) 



Therefore, the calculation of density matrix is largely de- 
termined by the evolution in the zero-quantum space. 
Under Hartmann-Hahn match condition uju = Wis, the 
Hamiltonian in the zero-quantum space becomes 



U A (t) = d(t)a 



A 

z ' 



(8) 



which is time dependent. For every particular time, it 
can be represented by a zero-order average Hamiltonian 

A ft 

n A( -°\t) = ^ d{t')dt' = d(t)a A 
t Jo 

= — {2V2sin(2f3)[sin(oj r t + 7) — 5*71(7)] ^ 

- sin 2 {/3)[sin{2uj r t + 2j) - sin(2~f)]}a A . 

Since [H A (tf), H A (t")} = for any t' and t" , all the 
higher order Hamiltonians vanish. Consequently, the 
zero-order average Hamiltonian represents exactly the 
Hamiltonian of the spin system. It turns out that the 
average Hamiltonian itself is time dependent, but for any 
particular time t it can be treated as a time independent 
Hamiltonian in the evolution of density matrix. There- 
fore, for a given average Hamiltonian the density matrix 
in the zero-quantum space can be derived 



a A (t) = e- id ^<<7 A e a W ff . 



cos(d(t))a A — sin(d(t))a A . 



(10) 



Base on the Eqs. ([6]), (0, and (|T0|). the evolution of 
density matrix becomes 

p{t) = a A (t)+^{t) 

— cos(d(t))a A — sin(d(t))a A + of, 

which leads to a CP dynamics 

CP{t) = Trace < S y ■ p(t) > 



(11) 



= Trace < (of - a A ) • p(t) > 
= 1 - cos(d(t)). 



(12) 



The CP dynamics of powder sample can be derived by 
integration over the Euler solid angle on the sphere 

{1 - cos[d(t)]}sin(f3)d(3d-f. (13) 

Considering the effects of I spin diffusion and spin lat- 
tice relaxation in the rotating frame, the CP dynamics 
becomes 



M(t) = M -{l-^e 
-fl.it 



-Ri 



(14) 



e- Hlt {l-2-CP(t))}-e T ip, 



where R is the spin diffusion rate of I spin, R\ is a rate 
which results in the oscillation damping by the remote I 
spins, and T\ p is the spin lattice relaxation in the rotating 
frame of the I spins. 

When oj r -> 0, the Eq. (|T4"|| reduce to: 



M{t)=M .{l-\e- Ri - 1 -e 



-flit 



^cosidt)} ■ e T ip (15) 



which, except for the Tl effect, is identical to the result 
by Miiller et al '3\ for a stationary sample. 

In the above calculation, the offsets, Auj and Awg for 
I and S spin, are ignored. Under off-resonance condi- 
tion, the Hartman-Hahn match condition is determined 
by effective fields, i.e. wu e — uj\Se, where 



and 



UJiSe = A/Awl 



(16) 



Correspondingly, the dipolar interaction is split up into 
two parts (perpendicular and parallel) 



d(t)I z S z = sin(0 I )sin(0 s )d(t)I ze S ze 
+ cos{9 I )cos(9s)d(t)I ye Sy e , 



where 



?/ = cos 1 (Awi/w 1I e) 

9 S = COS~ 1 (Au!s/u!lSe)- 



and 



(17) 



(18) 



In Eq. (|17l) . the perpendicular term is responsible for 
polarization transfer while the parallel term is a small 
perturbation which can be ignored here. In the perpen- 
dicular term, the dipolar interaction appears to be scaled 
by sin(9i)sin(8s)- It in turn cases a scaled polarization 
transfer rate by the same factor. In this case, the initial 
polarization should be placed along the effective filed of I 
spin, and the polarization of S spin is built up along the 
effective field of S spin. 



3. COMPUTER SIMULATION 

A home-made Java program termed QCNMR (quan- 
tum computation of NMR) , which is based on the evolu- 
tion of density matrix, is used for comparing the results 
from the analytical solution and computer simulation. 
All the solid lines shown in FIG H] and FIG 12] are derived 
from the analytic solution (Eq. (TT2"j) and (Q3|)) while the 
solid circles in the figures are given by the computer sim- 
ulation and NMR experiment, respectively. 

The CP dynamics with a MAS speed of 2 kHz is shown 
in FIGfTJ It appears to be periodic with a period which is 
the same as the period of MAS (T = 2ir/uj r ). The polar- 
ization of initial buildup and two nulls within the period 
(T = 500us) are caused by the interference between dipo- 
lar oscillation and MAS. This pattern is unique in slow 
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FIG. 1: Comparison of CP dynamics derived from the ana- 
lytical solution (solid line) and QCNMR (solid circles) with 
the conditions: ujh/2it = uis/27r=80kHz and heteronuclear 
dipolar constant d/-7r=5kHz. 
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4. CONCLUSION 

Under Hartman-Hahn match condition, CP dynam- 
ics can be derived analytically in the zero- and double- 
quantum spaces. The solution is valid for any MAS speed 
and offsets. In particular, the dynamics for a stationary 
sample appears when MAS speed approaches zero. As 
many other methods, the CP dynamics provides valu- 
able molecular structural information. Similar to REDO 
experiment [2l| . this analytic solution also provides a 
measure of dipolar coupling constant (or distance) for a 
strongly coupled system that is surrounded by a moder- 
ately coupled network. Unlike many other others 
[26j], this method does not required high MAS speed. If 
LG spin locking is applied to proton channel, the result 
of Least-Square fitting in FIG|2]will be better because all 
homonuclear coupling is decoupled. 

For a time dependent system, it is unlikely to find a 
systematic way for analytic solutions. This may explain 
why CP dynamics discussed in this article was delayed for 
so long a time. For an inhomogeneously broaden system 
([H(t),H(t")] = 0) [13, all high order average Hamil- 
tonians become zero except for the zero order average 
Hamiltonian that can be calculate conveniently. This 
method is quite general and can be used in NMR, op- 
tics, quantum computing and quantum mechanic related 
problems of a similar nature. 



FIG. 2: CP dynamics from a powder sample of alanine derived 
from the analytical solution (solid line) and experiments with 
the conditions: R~ x = 290.8fj.sec, R^ 1 — I37.9fisec, Tip = 
1.867msec, luh /2tt = k>i,g/27r=80kHz, and the distance be- 
tween 1 H and 13 C spin is 1.09A [la| . The experiment was 
done by Bruker Avance 300MHz NMR Instrument. 



MAS speed. In the above calculation R, R±, and T\ p 
in the equation are all neglected for better comparison 
between theory and simulation. 

In FIGlH we show the dynamics of a powder sample 
under a MAS speed of 5 kHz. Unlike a single crystal, 
the oscillation is strongly damped by the orientations of 
heteronuclear dipolar tensors. The polarization increases 
gradually as the CP and spin diffusion take place. The 
experiment results are normalized according to quantita- 
tive CP experiment with a reciprocity relation [l9|, , 
while the solid line is normalized by the Eq. (fl"4")) . In this 
case the match of two curves depends not only on the 
patterns but also on the specific values as well. 

It can be seen from FIGQ]and FIG|H au the simulated 
results agree well with the analytical solutions, demon- 
strating the validity of the analytical solution. The only 
approximation in the derivation is -\-ui\s\ S> 
which is well fulfilled in practice. 
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